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Abstract 



Precise numerical calculations of bound states of three-atomic He- 
lium cluster are performed. The modern techniques of solution of Fad- 
deev equations are combined to obtain an efficient numerical scheme. 
CO | Binding energies and other observables for ground and excited states 

are calculated. Geometric properties of the clusters are discussed. 



OS, 

1 Introduction 

• »h \ Small clusters of Helium attract the attention of specialists in different fields 

j>^' of physics. Fine experimental techniques are developed to observe these 

clusters [I], ^], [3]]. Different quantum chemistry approaches are used to pro 



duce numerous potential models of He-He interaction [17|, [18], [19], £(], [21], f22 
Model-free Monte-Carlo calculations were performed to check the accuracy 
of the models ||23||. The special attention is payed to three-body Helium 



clusters || ||, [7], §, || [10], [ll| because of their possible near-Efimov behavior 
[|]]. Complicated shape of the model potentials also makes Helium trimer a 
perfect touchstone for the computational methods of three-body bound state 
calculations | I3 |. 
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Although the investigation of He3 lasts already more than 20 years [||], 
some important physical questions has not received definite answer yet. One 
of the questions is dealt with speculations on Efimov-like nature of the He3 
bound states: how many excited states are supported by the best known 
model interactions? Can one estimate any differences in the number of bound 
states varying the model potentials being limited by the accuracy of contem- 
porary models? Another important question is dealt with the characteristics 
of He3 bound states. Can He3 trimers influence the results of experimental 
measurement of He2 dimer characteristics? To answer these questions one 
should know such important characteristics of the He3 cluster as the mean 
square radius of different states of the trimer and its geometric shape. 

In this paper we investigate the 4 He trimer performing direct calculations 
of 4 He3 bound states with different He-He model potentials. We base our 
calculations on Faddeev equations in configuration space because of simplic- 
ity of numerical approximation of Faddeev components in comparison with 
a wave function. In the case of Faddeev equations the boundary conditions 
are also much simpler. 

In the following sections the equations we have solved numerically are 
described, some observables for different states of He3 and for He2 are pre- 
sented. 



2 Faddeev equations for bound states 



According to the Faddeev formalism |L4] the wave function of three particles 



is expressed in terms of Faddeev components $ 

^(xi,yi) =$i(xi,yi) + $2(x 2 ,y 2 ) + $ 3 (x3,y 3 ) , (1) 

where x Q and y a are Jackobi coordinates corresponding to the fixed pair a 
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Here r a are the positions of the particles in the center-of-mass frame. The 
Faddeev components obey the set of three equations 

(-A x -A y + V a (x a ) - £)$ a (x a , y Q ) = -V a {x a ) Y, "Mx/3, Y/s) 

fca , (3) 

a = 1,2,3 

where V a (x. a ) stands for pairwise potential. To make this system of equa- 
tions suitable for numerical calculations one should take into account the 
symmetries of the physical system. Exploiting the identity of Helium atoms 
in the trimer one can reduce the equations ([3|) to one equation [|14[]. Since 



all the model potentials are central it is possible to factor out the degrees 



of freedom corresponding to the rotations of the whole cluster \T5\. For the 



case of zero total angular momentum the reduced Faddeev equation reads 

'' - — (- -) — (I :-)" — + 

(4) 
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C + and C~ are cyclic and anticyclic permutation operators acting on the 
coordinates x, y and z as follows 
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C±x = {^ + ^T^xy Z ) l '\ 
C±y=(^ + ^±fxyz)V\ 

p± ^4 + 4 2 X i)^ 

C^xC^y 
The asymptotic boundary condition for bound states consists of two terms 

x e -k 3 (x 2 +y 2 )l 

$(x,y,z)~ (f) 2 {x)e~ kyy + A(-,z) 
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where <fe(a;) is the two-body bound state wave function, k y = \JE 2 — E 3 , 
k 3 = v / — E 3 , E 2 is the energy of the two-body bound state and E 3 is the 
energy of the three-body system. The second term corresponding to virtual 
decay of three body bound state into three single particles decreases much 
faster than the first one which corresponds to virtual decay into a particle 
and two-body cluster. In our calculations we neglect the second term in the 
asymptotics introducing the following approximate boundary conditions for 
the Faddeev component at sufficiently large distances R x and R y 

d x ${x,y,z)l x=Rx _ r- 

K2 = *y 2 ' 



d y ^x,y,z)[ y=Ry 
<5>(x,y,z)[ y=Ry 



To calculate the bound state energy and the corresponding Faddeev com- 
ponent one has to solve the equation (f|) with the approximate boundary 
condition ([|) . The numerical scheme we have chosen to perform the calcula- 
tions is based on tensor-trick algorithm |16|] . In this paper we do not describe 



the realization of the numerical methods exploited but only underline some 
essential features of our approach. They are 



1. total angular momentum representation fL5| , 

2. tensor-trick algorithm [IE], 



3. Cartesian coordinates [[K]]. 

The total angular momentum representation itself is a strong method of par- 
tial analysis allowing to take into account contribution of all the angular 



momentum states of two-body subsystems at once [pl|. Tensor-trick algo- 



rithm [Tj| is known to be a powerful method of solution of Faddeev equations 
for bound states. Being applied to the equations in total angular momentum 
representation it leads to effective computational scheme which makes pos- 
sible to use all the advantages of Cartesian coordinates. In particular using 
Cartesian coordinates [TIJ one can obtain a criterion to select the optimal 



grid in the coordinate x. This criterion comes from the asymptotic behavior 
of Faddeev component 



$(x,y,z) ~ <p2(x)e 



-k y y 



where <fi2{x) is the two-body bound state wave function. That is why the 
properly chosen grid in x should support the correct two-body wave func- 
tion. Comparing the binding energy of two-body subsystem calculated on the 
"three-body" grid with the exact results one can estimate the lower bound 
for a numerical error of a three-body calculation. 

Thus the usage of total angular momentum representation has allowed 
us to construct an efficient numerical scheme combining the most advanced 
methods proposed during the last decade. 

3 Results of calculations 

Having the equation (|4]) solved numerically one has the value of the en- 
ergy of 3-body state E 3 and the corresponding Faddeev component $>(x,y,z) 
for a particular model potential. Comparing the observables calculated for 
different potential models one can estimate the bounds limiting the values 
of these observables for the real system. Eight different potential models 
were used in our calculations: HFD-HE2 pf, HFD-B(He) p|, LM2M1 @ 



LM2M2 0, HFDID gg, LM2M1 and LM2M2 without add-on correction 
term p0 |, TTYPT [22 1 . In the Tab. 1 we give the values of trimer energies 
for ground and excited states. To confirm the accuracy of our calculation 
we also present the values of dimer binding energy calculated on the grid 
used in three-body calculations Ei and the exact results E%. The difference 
between these values can be regarded as the lower bound for the error of 
our approach. In the Tabs. 2 and 3 we demonstrate the convergence of the 
calculated energies with respect to the number of grid points used in the 
calculations. The results of other authors for the most known potentials are 
given in the Tab. 7. The best agreement is observed with the results of [BJ 



and |] T| . In the ref . || no angular momentum cut-off is made, that makes it 
the closest one to our approach. In all other papers some kind of partial wave 
decomposition is performed and finite number of angular basic functions is 



taken into account. The most complete basis is used in the ref. |]ll|. The 



agreement between our calculations and the result of |TTJ for the excited state 



is impressive, but the ground state energy of |L1]] is about one percent less 
than our result. Consideration of the geometric properties of Helium trimer 
can clarify the possible nature of this difference. 

Since the Faddeev component is calculated, the wave function can be 



recovered as follows 



ip{x, y, z) = $(x, y, z) + xy ( — — h 



x^y^ x y 

where x ± = C ± x and y^ 1 = C ± y. Having the wave function recovered one 
can investigate the shape properties of the system. The most intuitive way 
to visualize the results of the calculations is to draw a one-particle density 
function defined as 

p(ri) = / dr 2 dr 3 |^(ri,r 2 ,r 3 )| 2 , 

where 

il){x(r 1 ,r 2 , r 3 ), y(r x , r 2 , r 3 ), z(r 1 ,r 2 , r 3 )) 



*(ri,r a ,r 8 ) 



47r 3 x(ri, r 2 , r 3 )y(ri, r 2 , r 3 ) 



the functions x(ri, r 2 , r 3 ), y(r l5 r 2 , r 3 ) and 2(1*1, r 2 , r 3 ) are defined according 
to (0) and (|), the function ip(x,y,z) is normalized to unit. Due to the 
symmetry of the system the one-particle density function is a function of the 
7*1 = I i"i I coordinate only. Taking into account the relation ri = Voyi we 
get 

p(r) = Vo / dxdz \tp(x, —=, z)\ 2 . 

Omitting the integration over z we define a conditional density function 
p(r, z) that presents a spatial distribution for the particle 1 when the other 
two particles are located along the fixed axis. It is useful to plot this func- 
tion in coordinates (ri, r a ) such that r; = rz is a projection of the particle 1 
position to the axis connecting the other particles and r a = t§|T(1 — z 2 )^ is 
a projection to the orthogonal axis. Three-dimensional plots of the function 
p((rf + r 2 ) 1 / 2 , cosarctan— ) corresponding to the ground and excited states 
of the trimer calculated with LM2M2 potential are presented on the Fig. 1 
and Fig. 2. The conditional density function of the ground state decreases 
democratically in all the directions. The density function of the excited state 
has two distinguishable maximums demonstrating the linear structure of the 
cluster. This structure has a simple physical explanation. The most probable 
positions of a particle in the excited state lie around two other particles. At 
the distances where two particles are well separated the third one forms a 
dimer-like bound state with each of them. This interpretation agrees with 
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the clusterisation coefficients presented in the Tab. 4. These coefficients are 
calculated as a norm of the function f c defined as follows 



fc(y) = / dxdz$(x,y,z)(j)2(x) , 



where 4>2{x) is the dimer wave function. The values of ||/ c (?/)|| 2 shown in the 
Tab. 5 demonstrate the dominating role of a two-body contribution to the 
trimer excited state whereas in the ground state this contribution is rather 
small. We could suppose that this dominating contribution of the cluster 
wave in the excited state has ensured fast convergence of the hyperspherical 



adiabatic expansion in the paper [II] to the correct value, but to get the 



same order of accuracy for the ground state possibly more basic functions 
should be taken into account. 

Very demonstrative example of the advantage of Faddeev equations over 
the Schroedinger one in bound-state calculations is given in the Tabs. 8 and 9. 
Here we present the contribution of different angular states to the Faddeev 
component and to the wave function calculated as 

C n =\\f n (x,y)\\ 2 
f n (x,y)=j dzF(x,y,z)P n (z) , 

where P n (z) are the normalized Legendre polynomials, F(x, y, z) is the Fad- 
deev component or the wave function, n = 0,2,4. The angular coefficients 
for the Faddeev component decrease much faster than the wave function co- 
efficients. The Tab. 8 also demonstrates that more angular functions should 
be taken into account in the ground state calculations. 

4 Conclusions 

The high accuracy calculations of He3 bound states were performed on the 
base of the most advanced few-body calculations techniques. Eight different 
potential models were used. For every potential model, either more (LM2M2, 
TTYPT) or less realistic one (LM2M2a, HFD-ID), two and only two bound 
states are found. The properties of these states are very different. The 
ground state is strongly bound, whereas the binding energy of the excited 
state is comparable with the binding energy of dimer. The sizes of these two 
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states also differs much. The characteristic size of the ground state either 
estimated by (r) or (r 2 ) 1 / 2 (Tabs. 5 and 6) is approximately 10 times less 
than the size of dimer molecule, but the size of the excited state has the 
same order of magnitude with the dimer's one. This estimation shows the 
necessity to check for the absence of trimers in the experimental media during 
the measurement of dimer properties and vice versa. 
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Table 1: The energy of the He2 and He3 bound states 



Potential 


E 2 ,mK 


E 2 , mK 


E 3 , K 


E*, mK 


HFD-A 


-0.830124 


-0.8305 


0.11713 


1.665 


HFD-B 


-1.685419 


-1.68540 


0.13298 


2.734 


HFD-ID 


-0.40229 


-0.4024 


0.10612 


1.06 


LM2M1 


-1.20909 


-1.212 


0.12465 


2.155 


LM2M2 


-1.303482 


-1.304 


0.12641 


2.271 


LM2Mla 


-1.52590 


-1.527 


0.13024 


2.543 


LM2M2a 


-1.798436 


-1.795 


0.13471 


2.868 


TTY 


-1.312262 


-1.3121 


0.12640 


2.280 



Table 2: Convergence of the He3 excited state energy with respect to the 
number of gridpoints 



Grid 


E*,A~ 2 x 10~ 5 


E 2 ,A- 2 x 10~ 5 


45x45x9 


-22.819 


-14.123 


60x60x9 


-22.568 


-13.913 


60x60x15 


-22.570 


-13.913 


75x75x9 


-22.561 


-13.907 


75x75x15 


-22.563 


-13.907 


90x75x9 


-22.567 


-13.912 


105x75x9 


22.555 


-13.902 



Table 3: Convergence of the He3 ground state energy with respect to the 
number of gridpoints 



Grid 


E 3 ,A~ 2 x 10~ 5 


E 2 ,A~ 2 x 10~ 5 


45x45x15 


-1096.35 


-13.839 


60x60x9 


-1096.72 


-13.894 


60x60x15 


-1097.11 


-13.894 


60x60x21 


-1097.11 


-13.894 


105x60x18 


-1097.19 


-13.9062 


105x75x15 


-1097.25 


-13.9062 
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Table 4: Contribution of cluster wave to the Faddeev component 



Potential 


ll/e|| 2 


H/cT 


HFD-A 


0.2094 


0.9077 


HFD-B 


0.2717 


0.9432 


HFD-ID 


0.1555 


0.8537 


LM2M1 


0.2412 


0.9283 


LM2M2 


0.2479 


0.9319 


LM2Mla 


0.2624 


0.9390 


LM2M2a 


0.2780 


0.9458 


TTY 


0.2487 


0.9323 



Table 5: The mean square radius of Helium molecules, A 



Potential 


Ground state of He3 


Excited state of He3 


He 2 


HFD-A 


6.46 


66.25 


88.18 


HFD-B 


6.23 


57.89 


62.71 


HFD-ID 


6.64 


75.38 


126.73 


LM2M1 


6.35 


61.74 


73.54 


LM2M2 


6.32 


60.85 


70.93 


LM2Mla 


6.27 


59.03 


65.76 


LM2M2a 


6.21 


57.17 


60.79 


TTYPT 


6.33 


60.81 


70.70 
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Table 6: The mean radius of Helium molecules, A 



Potential 


Ground state of He3 


Excited state of He3 


He 2 


HFD-A 


5.65 


55.26 


64.21 


HFD-B 


5.48 


48.33 


46.18 


HFD-ID 


5.80 


62.75 


91.50 


LM2M1 


5.57 


51.53 


53.85 


LM2M2 


5.55 


50.79 


52.00 


LM2Mla 


5.51 


49.28 


48.34 


LM2M2a 


5.46 


47.72 


44.82 


TTYPT 


5.55 


50.76 


51.84 



Table 7: Comparison with the results of other authors 



Observable 


This work 


1 


11 


11 


|10| 


HFD-A(He) 


E 3 , K 


-0.1171 


-0.117 


-0.114 


-0.11 


-0.107 


E*, niK 


-1.665 




-1.74 


-1.6 




HFD-B(He) 


E 3l K 


-0.1330 




-0.132 




-0.130 


E*, mK 


-2.734 




-2.83 







LM2M2 


Observable 


This work 


|H] 


|12| 


E 3 , K 


-0.1264 


-0.1252 


-0.219 


E*, mK 


-2.271 


-2.269 


-1.73 


< r 2 > 1 /2, A 


6.32 


6.24 


7.4 


< rl > 1 /2, A 


60.85 


60.86 


50.3 
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Table 8: Contribution of different two-body angular states to the Faddeev 
component 





Ground state 


Excited state 


Potential 


S 


D 


G 


S 


D 


G 


HFD-A 


0.9991043 


0.0008859 


0.0000095 


0.9999964 


0.0000035 


0.0000000 


HFD-B 


0.9990000 


0.0009890 


0.0000107 


0.9999952 


0.0000048 


0.0000001 


HFD-ID 


0.9991709 


0.0008200 


0.0000088 


0.9999972 


0.0000028 


0.0000000 


LM2M1 


0.9990505 


0.0009390 


0.0000101 


0.9999958 


0.0000042 


0.0000000 


LM2M2 


0.9990393 


0.0009500 


0.0000103 


0.9999957 


0.0000043 


0.0000000 


LM2Mla 


0.9990129 


0.0009762 


0.0000105 


0.9999954 


0.0000046 


0.0000001 


LM2M2a 


0.9989834 


0.0010053 


0.0000109 


0.9999950 


0.0000049 


0.0000001 


TTY 


0.9990332 


0.0009561 


0.0000104 


0.9999956 


0.0000043 


0.0000000 



Table 9: Contribution of different two-body angular states to the wave func- 
tion 





Ground state 


Excited state 


Potential 


S 


D 


G 


S 


D 


G 


HFD-A 


0.95416 


0.03198 


0.00877 


0.90957 


0.07543 


0.01331 


HFD-B 


0.95193 


0.03365 


0.00947 


0.89710 


0.08546 


0.01546 


HFD-ID 


0.95493 


0.03116 


0.00905 


0.91919 


0.06763 


0.01170 


LM2M1 


0.95323 


0.03277 


0.00891 


0.90337 


0.08043 


0.01437 


LM2M2 


0.95303 


0.03294 


0.00893 


0.90201 


0.08152 


0.01460 


LM2Mla 


0.95259 


0.03332 


0.00899 


0.89904 


0.08391 


0.01512 


LM2M2a 


0.95210 


0.03374 


0.00906 


0.89574 


0.08654 


0.01569 


TTY 


0.95245 


0.03318 


0.00941 


0.90186 


0.08164 


0.01463 
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Figure 1: Conditional one-particle density function of the He3 ground state, 
r h r a in A 
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Figure 2: Conditional one-particle density function of the He 3 excited state, 

r h r a in A 




Figure 3: He 3 ground state density function, r in A 
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Figure 4: He3 excited state density function, r in A 
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Figure 5: He2 dimer density function, r in A 
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